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One problem facing the Aerospace Environmental Support Center (Det 1, 4th 
WW, AWS) is the prediction of ionospheric effects related to energetic proton 
emissions from the sun. The intensity and duration of these charged particle 
emissions are related to the optical, radio, and X-ray emissions from the solar 
flares. Radio and X-ray data give an important indication of proton produc- 
tion. Cm radio bursts with a high flux density and a slope reflecting greater 
flux values at higher frequencies are often proton productive. Examination of 
the data summaries has shown that the meter classification of flares as type IX 
or type IV, even though qualitative, should be used in conjunction with the cm 
radio criteria. The peak intensity o of the proton flux E >_ 10 meV correlates 
better with the integrated 1 to 8 A X-ray bursts than with the integrated cm 
radio bursts, peak cm radio flux, or peak X-ray flux. Free-free transition 
theory was used to relate the X-ray flux. Free-free transition theory was used 
to relate the X-ray intensities recorded by different satellites with the peak 
proton fluxes observed at the earth. 

Both the anisotropic and isotropic diffusion theories could be used to ex- 
trapolate proton fluxes for E ^ 10 meV for over 50% of the particle events. 
The isotropic diffusion theory uses a diffusion coefficient. 

D - Mr 0 

It was found that M and 8 tended to be functions of flare position on the 
solar disk. 

A measurement of the interplanetary flux in near earth space gives a good 
indication of the polar cap fluxes. It was found that the 30 MHz absorption 
over the poles during a PCA is proportional to the square root of the integral 
proton flux E il me v i n interplanetary space. 

J * KA I. 2 

with K * 8 + 2 and J in protons/cm^-sec-ster . 

Without the moral support of the personnel assigned to the Solar Forecast 
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ten. My thanks go to Mr. Ray Cormier of Air Force Cambridge Research Labora- 
tories for allowing me to use his riometer data. The Vela X-ray data were 
kindly supplied by Dr. Kunz, Dr. Conner, and Dr. Bame of Los Alamos Scientific 
Laboratories. The author is solely responsible for the opinions and material 
presented in this report. The consent of the above-mentioned scientists to use 
their data does not necessarily imply their agreement with the results. 


I . Introduction 

One problem facing the Air Weather 
Service Solar Forecast Center is the pre- 
diction of ionospheric effects related to 
energetic proton emissions from the sun. 
The intensity, and duration of these 
charged particle emissions are related to 
the optical, radio, and X-ray emissions 
from the solar flares. Data summaries 


appearing in pre-published form, such as 
the Explorer 34 proton results in the 
Solar-Geophysical Data bulletins published 
by ESSA (ref. 1) and the Geophysics and 
Space Data ’ Bulletin by AFCRL (ref. 277 
have allowed the evaluation of many impor- 
tant parameters . 

The work reported in the paper was 
done for an operational unit. Therefore, 


356 
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the data base available to the operational 
solar forecaster was used. This was not 
the most precise data base available. 
Several simplifying assumptions had to be 
made which will become apparent in this 
paper. Nevertheless, the results clearly 
indicate that short-term proton event pre- 
dictions are possible and indicate certain 
limitations on several current theories. 


j « N P 

This paper presents some experimental de- 
terminations of N and ?. The actual pre- 
diction of proton production has been 
previously published (ref. 3). 

III. Proton Acceleration and X-Radiation 
Flux 


An interdisciplinary approach to PCA 
prediction and ionospheric forecasting 
allows some effects to be calculated from 
first principles shortly after the optical 
flare occurs. Solar radio and X-ray burst 
observations are useful in determining if 
protons were actually produced during the 
flare. These bursts also give an indica- 
tion of the number of protons produced. 
The position of the flare on the solar 
disk and the solar wind velocity can be 
used to indicate the risetime and the peak 
intensity of the event in the near earth 
space. Using the first few proton data 
points, the proton intensities can be ex- 
trapolated into the future. Ground-based 
magnetic measurements can then be used to 
indicate the spatial extent of the precip- 
itation and its uniformity. Satellite 
data and real time magnetic measurements 
can be used to update the forecasts during 
the event. On the basis of the extrapola- 
ted fluxes, the next step is to calculate 
the ionization caused by the protons using 
an atmospheric model and a two-ion D- 
region model for PCA events to obtain the 
electron profiles. Finally, the absorp- 
tion on vertical and oblique paths can be 
calculated for high-frequency communica- 
tion systems across the polar regions. 
Real time satellite and riometer measure- 
ments can be used to verify the prediction. 

This paper covers several of the dif- 
ferent phases of PCA prediction. First, 
the X-ray characteristics of proton pro- 
ducing flares are examined. Then the two 
theories relating to particle propagation 
in the interplanetary medium are investi- 
gated for their applicability. Finally, 
the relationships between 30 MHz riometer 
absorption and the exospheric proton flux 
are discussed. The problem of electron 
profiles and oblique absorption will not 
be covered in this report. The radio 
characteristics of the flare and the spa- 
tial extent of the polar cap absorption 
event will not be covered because they 
have been discussed previously (ref. 3) . 

II. General Discussion 

In order to meaningfully predict pro- 
ton fluxes at the earth before the arrival 
of large numbers of protons, one must use 
the electromagnetic characteristics of the 
parent flare. The simplest set of assump- 
tions one can make about the peak proton 
fluxes are that (1) the peak flux is pro- 
portional to the number of particles 
accelerated by the parent flare, N, and 
(2) the transport of the particles to the 
earth is influenced by the interplanetary 
medium, P. Thus, the flux, J is just 


There are several methods of deter- 
mining N , such as integrated or peak radio 
intensity, H-a flare classification, and 
integrated or peak X-ray intensity. Of 
these, the integrated X-ray intensity 
seems to be the best quantitative measure 
of proton acceleration. If the number of 
electrons accelerated during the flare is 
proportional to the number of protons ac- 
celerated, and if the intensity of X rays 
in some wavelength band is proportional to 
the number of electrons being accelerated, 
then the integrated intensity of X-radia- 
tion in that band ought to be proportional 
to the number of protons accelerated. For 
simplicity, it was assumed that the X-ray 
bursts had exponential rise and exponen- 
tial decay time. Then 
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where t r and t^ are the exponential rise 
and decay times , I ma x is the maximum X-ray 
intensity, and I± n t is the integrated X- 
ray intensity. The use of such a simple 
relationship allows a forecaster to make 
predictions with minimum information. 

The largest, most consistent X-ray 
data base available is the Vela X-ray data 
(refs. 5, 6). The 0.5 to 5 A X-ray data 
were integrated using equation (3) . The 
result of plotting the integrated X-ray 
intensity versus the peak 10-meV proton 
fluxes observed by IMP F (Explorer 33) is 
given in figure 1 (ref. 1). Even though 
these points are not corrected for propa- 
gation through the interplanetary space, 
the correlation is striking. 


One must generalize these results 
shown in figure 1 , to those obtainable 
from different sensors. This was a prob- 
lem since the integrated X-ray flux pre- 
sented in figure 1 was a measure of the 
flux deposited inside the detector after 
transmission through detector windows. 
The most commonly used assumption when de- 
termining the flux incident on an X-ray 
detector is that the photon spectrum is a 
blackbody Planckian spectrum of the form 


W(X ,T)dX 


87rhc dX 

X 5 e hc/XkT -l 
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Figure 1. Integrated X-ray flux deposited 
in 0.5 to 5.0 A Vela launch 2, 3, 4 de- 
tectors vs. peak proton flux (E >_ 10 
meV protons/cm^-sec-ster) . Propagation 
factors have not been removed. 


In order to compare the data obtained 
from different satellites, this spectral 
form was folded into the different detec- 
tor responses as a function of solar flare 
temperature. When this form was used, 
different detectors did not give consis- 
tent results. Figure 2 shows the flare 



Figure 2. Solar flare temperature Vela 
launch 2, 3, 4 vs. solar falre tempera- 
ture Vela launch 5. Figure 2a is the 
grey body assumption while figure 2b is 
the free-free transition assumption. 


temperatures determined by two different 
methods from two distinctly different sets 
of satellite instrumentation. Note that 
the grey body assumption (fig. '2a) does 
not fit the data accurately. Furthermore, 
when the grey body assumption was employed 
the correlation between the 1 to 8 A inte- 
grated X-ray flux and the peak proton flux 
disappeared . 

The functional form of the X-ray en- 
ergy spectrum which allowed the correla- 
tion between the 1 to 8 A X-ray flux and 
the peak proton flux to remain and which 
better correlated the solar flare tempera- 
tures given in figure 2b was 

” U ' T > * ^ -- - ric/XkT 
a e 

This is the functional form obtained by 
assuming the X rays are produced by a Max- 
wellian distribution of electrons undergo- 
ing free-free transitions (i.e., brems- 
strahlung ) (ref. 7, 8). This assumption 
allowed the peak proton flux to be corre- 
lated with the 1 to 8 A X-ray flux, which 
is given in figure 3. Not enough data 
have been published at the present time to 



Figure 3 . Integrated X-ray flux incident 
on a detector in the 1 to 8 A band vs. 
proton flux E > IQ meV. Points were 
derived from figure 1 using the free- 
free transition assumption. Triangle 
points were derived assuming the aver- 
age proton flare temperature of 7xio 6 °K. 


allow the complete verification of this 
analysis. However, the few points which 
have been obtained from Explorer 44 are 
presented in figure 1, but the scatter is 
well within that expected. The average 
solar flare temperature for the ' proton 
producing flares was 7><10 6 °K. This tem- 
perature was obtained from the Vela launch 
2, 3, 4 sensors using free-free transition 
theory . 
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IV, Arrival of Protons 


D = Mr g 


(4) 


There are two theoretical frameworks 
within which the proton observations could 
be analyzed. The first is the anisotropic 
diffusion with boundary theoretical frame- 
work (ADB) of Burlaga (ref, 9), while the 
second is the isotropic diffusion theoret- 
ical framework (refs, 10, 11). Both theo- 
ries have definite shortcomings. These 
shortcomings are illustrated in figure 4. 
If the theoretical models were adequate at 
10 meV, one would expect systematic trends 
when plotting the proton risetimes against 
the 1/e decay times. The scatter in this 
diagram illustrates the problems involved 
in attempting to fit the experimental 
data. Similar plots for E > 30 and E > 60 
meV protons do show more systematic re- 
sults . 



Figure 4. The time in hours it took 10- 
meV protons to peak as a function of 
the time it took for the protons to 
reach 1/e of their maximum value. 


The anisotropic diffusion with bound- 
ary (ADB) theory assumes that the inter- 
planetary magnetic field lines between the 
sun and the earth are spirals with irregu- 
larities which can effectively scatter the 
solar protons. The solar proton flux in 
this theory will be anisotropic and the 
flux will be greatest in the direction of 
the field. Beyond 1 A.U. , the scattering 
centers change so that protons cannot be 
scattered back into the diffusing region 
once they have reached the transition re- 
gion. Thus, in this model, it is assumed 
that there is an absorbing boundary at 
some heliocentric distance greater than 
1 A.U. Since systematic trends did not 
appear in applying this theory, and since 
this theory has been discussed in detail 
(ref. 3) , it will not be examined in any 
more detail in this paper. 

One can also analyze each proton 
event in terms of the isotropic diffusion 
theory (refs. 10, 11). In this theory, 
the diffusion coefficient, D, is given by 


where r is the heliocentric radial dis- 
tance in A.U. , and M and 3 are parameters 
which may and often do depend upon the 
particle energy E. The units of Mr^ " 2 
are (hours) . Krimigis shows that for 
E < 50 meV, the parameter M does depend 
upon the proton energy. Thus, if one 
wants the time behavior for 10 meV pro- 
tons, one must use a data base of 10 meV. 

It can be shown that 


t 


^ JL r 2 ^ 
max 3M (2-3) 


(5) 


where M and 3 are the coefficients given 
in equation (2) and t is the time it 
takes for the protons to x reach maximum in- 
tensity. Since t does depend upon the 
solar flare heliogra^hic longitude, it is 
possible to empirically determine M and 3 
as a function of solar flare position even 
though such an assumption violates the 
basic theoretical considerations. 

Approximately 80% of the observed 
events could be analyzed within the iso- 
tropic diffusion framework. Figure 5 
shows one of these events. This figure 
also shows that there is an area of agree- 
ment between the ADB theory and the iso- 
tropic diffusion theory. A certain amount 
of judgment was necessary in determining 6 
from the graphs of 1/t vs. In lt( 3 / 2 ~^' 
because of poor time resolution and uncer- 
tainties in count rates and particle back- 
grounds. This event also illustrates what 
happens when a proton event has considera- 
ble structure. Many events showed the 
characteristic hook at small values of 
1/t. The value of 3 used was that value 
which best fit the first several values of 
1/t and which also adequately fit the long 
time behavior at small values of 1/t. 



Figure 5. Intensity-time profile during 
26 April 1968 event for protons E ^ 10 
meV measured by Explorer 34. This 
event had considerable structure. 
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Figure 6 gives scatter plots of 3 and 
M as a function of 0 for 10 and 30 meV 
protons as a function of flare position. 
The 0 was determined in the ADB theoreti- 
cal framework (ref. 9). Several features 
do stand out. The low values of 3 and M 
for 0 less than 0.25 radian are probably 
due to the diffusion being one or two 
dimensional instead of three dimensional 
as assumed in this analysis. There does 
seem to be a trend as a function of 0 for 
0 > 0.5 radian. Although there is a large 
amount of scatter, the higher energy par- 
ticles tend to have larger values of M and 
3 than the lower energy particles. 
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Figure 6. The diffusion coef f icients r M 
and 3 as a function of flare position 
and the Earth-Sun Archimedes spiral 
angle 0 . The higher energy particles 
tend to have higher values of M and 3 . 


These results can be refined. Better 
time resolution will lead to less scatter 
in the data. The better the time resolu- 
tion, the easier it was to determine M and 
3 . Data for larger values of 0 have not 
yet become available with the necessary 
count rate accuracy or time resolution. 

V. Riometer Absorption 

Once one has predicted the intensity 
of protons in interplanetary space, one 
can predict an upper limit to the vertical 
30 MHz riometer absorption. 

On the basis of theoretical calcula- 
tions (ref. 12) , it can be concluded that 
the 30 MHz riometer absorption is propor- 
tional to the square root of the incident 
integral particle flux for E >_ 11 meV (ref. 


13) . An analysis of the Explorer 34 and 
Thule riometer data has helped confirm 
this result, 


J = KA 2 

where J is the flux in protons/cm -sec- 
ster, and A is the absorption in db of the 
30 MHz riometer. 

The Explorer 34 data and Thule riome- 
ter data have been examined. Figure 7 is 
a plot of the results for 10-meV particles 
One notes that these 10-meV particles 
scatter about the line of slope 0.5. The 
points tend to fall below the line at high 
absorption values and above the line at 
low absorption values. When the E ^ 20 
meV flux is plotted as the abscissa, the 
points at low absorption tend to fall be- 
low the line while those at high values 
tend to fall above the line. Thus, the 
square root of the absorption is directly 
proportional to the integral flux some- 
where between 10 and 20 meV. Within the 
limitations of this analysis, therefore, 
the integral flux for E >_ 11 meV is an ex- 
perimentally verified conclusion. The 
constant of proportionality, K, from this 
curve is 


K = 8 + 2 
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Figure 7. Riometer absorption at Thule as, 
a function of the interplanetary proton 
flux (E > 10 meV protons/cm 2 -sec-ster ) . 
If the square of the riometer absorp- 
tion were proportional to the flux, the 
points would fall on the straight line. 
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VI . Summary 

An interdisciplinary approach to pro- 
ton' event prediction allows some effects 
to be calculated from first principles 
shortly after the optical flare occurs. 
Solar X-ray burst characteristics are use- 
ful in determining if protons were pro- 
duced in sufficient quantity to give riom- 
eter absorption. The position of the 
flare on the solar disk and the solar wind 
velocity can be used to indicate the rise- 
time and peak intensity of the event in 
near earth space. Using the first few 
proton data points , the proton intensities 


can be extrapolated into the future using 
previously derived values for the proton 
diffusion coefficients given in figure 6. 
Satellite data can be used to update the 
predictions during the event by extrapola- 
ting the fluxes as a function of time. 
Experimental measurements have verified 
some of the assumptions and theoretical 
predictions for 30 MHz riometer absorption.. 
Thus, the satellite measurements can be 
employed to specify and predict proton 
fluxes and radio absorption. The approach 
taken in this report can be continued to 
give better predictions once more data be- 
come available. 
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